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We study the interaction between the Moon and the solar wind using a three-dimensional hybrid plasma solver. 
The proton fluxes and electromagnetical fields are presented for typical solar wind conditions with different 
magnetic field directions. We find two different wake structures for an interplanetary magnetic field that is 
perpendicular to the solar wind flow, and for one that is parallell to the flow. The wake for intermediate magnetic 
field directions will be a mix of these two extreme conditions. Several features are consistent with a fluid 
interaction, e.g., the presence of a rarefaction cone, and an increased magnetic field in the wake. There are 
however several kinetic features of the interaction. We find kinks in the magnetic field at the wake boundary. 
There are also density and magnetic field variations in the far wake, maybe from an ion beam instability related 
to the wake refill. The results are compared to observations by the WIND spacecraft during a wake crossing. The 
model magnetic field and ion velocities are in agreement with the measurements. The density and the electron 
temperature in the central wake are not as well captured by the model, probably from the lack of electron physics 
in the hybrid model. 
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1. Introduction 

Bodies that lack a significant atmosphere and internal 
magnetic fields, such as the Moon and asteroids, can to a 
first approximation be considered passive absorbers of the 
solar wind. The solar wind ions and electrons directly im¬ 
pact the surface of these bodies due to the lack of atmo¬ 
sphere, and the interplanetary magnetic field passes through 
the obstacle relatively undisturbed because the bodies are 
assumed to be non-conductive. Since the solar wind is ab¬ 
sorbed by the body, a wake is created behind the object. 
This wake is gradually filled by solar wind plasma down¬ 
stream of the body, through thermal expansion and the re¬ 
sulting ambipolar electric field, along the magnetic field 
lines. For a review of the Moon-solar wind interaction, 
including recent findings, see Halekas et al. (2011), who 
also point out the fact that we can view the near lunar space 
as a plasma laboratory. Although the Moon lacks any sig¬ 
nificant atmosphere and ionosphere there are many inter¬ 
esting features to this seemingly simple solar wind interac¬ 
tion. To understand the global interaction between the solar 
wind and the moon we need computer models that can give 
us a three-dimensional picture that complements the in-situ 
plasma observations. 

Magnetohydrodynamic (MHD) modeling of the Moon- 
solar wind interaction has a long history, starting with the 
work of Spreiter et al. (1970), who analytically treated 
the special case of an interplanetary magnetic field (IMF) 
aligned with the solar wind flow direction by finding a trans- 
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formation of the MHD equations into a hydrodynamic prob¬ 
lem. However, there are many kinetic processes that can¬ 
not be described by fluid models, e.g., the non-Maxwellian 
particle populations in the wake region, so particle mod¬ 
els should capture more of the physical processes. An 
approximation of the refill of the lunar wake is the gen¬ 
eral one-dimensional problem where plasma expands into 
a vacuum (Widner et al., 1971; Samir et al., 1983; Mora, 
2003), and such models have specifically been applied to 
the refill of the lunar wake (Farrell et al., 1998; Birch and 
Chapman, 2001). A nice property of such models is that 
even if they are one-dimensional, they can to some degree 
approximate a two-dimensional model of the lunar wake, 
where time corresponds to distance downstream the wake, 
i.e., the one-dimensional model is applied perpendicular to 
the solar wind flow and convects with the flow. A disad¬ 
vantage of one-dimensional models, e.g., along x, is that 
the magnetic field component along the x-axis is forced to 
be constant, from the requirement that V • B = 0. This 
is a severe restriction for modeling the lunar wake, where 
the field component across the wake can increase. Due to 
the computational complexity, particle in cell (PIC) models, 
that include electrons as particles, are for the foreseeable fu¬ 
ture limited to at most two spatial dimensions, so what has 
been investigated using such models is the interaction of an 
infinite cylinder with the solar wind (Harnett and Winglee, 
2003; Kimura and Nakagawa, 2008). 

Hybrid models, where ions are particles and electrons a 
fluid, are less computational demanding that full PIC mod¬ 
els and are suitable for studying processes of size compa¬ 
rable to the ion inertial length and gyro radius. A two- 
dimensional hybrid model by Travmcek et al. (2005) fo¬ 
cused mostly on wave activity in the lunar wake. The 
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Moon-solar wind interaction is however a fully three- 
dimensional problem, since the IMF component perpendic¬ 
ular to the solar wind flow introduces asymmetry into an 
otherwise cylindrical symmetric problem. A three dimen¬ 
sional hybrid model by Kallio (2005) had a fairly coarse 
grid on a small region, and did not handle the wake re¬ 
gion in a self consistent way (the electric field was specified 
in that region). Recently another three dimensional hybrid 
model has been presented in Wiehle et al. (2011). They 
compare observations by one of the ARTEMIS spacecrafts 
to a model run with time dependent inflow conditions, and 
find a good agreement. The time dependence enables com¬ 
parison with observations even during changing solar wind 
conditions, but makes it more difficult to see how the wake 
forms differently depending on the solar wind conditions. 

Here we also apply a three dimensional hybrid model 
to the Moon-solar wind interaction, but for steady inflow 
conditions. This allow us to study how the wake depends 
on the direction of the IMF, and in what follows we describe 
the method and present a comprehensive global overview of 
the different plasma quantities in the wake. 


2. Model 

In the hybrid approximation, ions are treated as particles, 
and electrons as a massless fluid. In what follows we use SI 
units. The trajectory of an ion, r(Q and v(Q, with charge q 
and mass m , is computed from the Lorentz force, 


dr d\ 

dt ’ dt 


— (E + v x B), 
m 


where E = E(r, t) is the electric field, and B = B(r, t) is 
the magnetic field. The electric field is given by 

E=-‘-(_j lxB + M - 1 (VxB)xB-V/ ?e ), (1) 

Pi 


where p\ is the ion charge density, Ji is the ion current, p Q 
is the electron pressure, and /jlq = An • 10 -7 [H m -1 ] is 
the magnetic constant. We assume that the electrons are 
an ideal gas, then p & = n Q kT Q , so the pressure is directly 
related to temperature (k is Boltzmann’s constant). There 
are several ways to handle the electron pressure (p. 8790, 
Winske and Quest, 1986). Here we assume that p Q is adi¬ 
abatic (small collision frequency), with an adiabatic in¬ 
dex, y = 5/3. Then the relative change in electron pres¬ 
sure is related to the relative change in electron density by 
P& OC |p e | K (Winske and Quest, 1986), and we have that 



where the zero subscript denote reference values (here the 
solar wind values at the inflow boundary). From charge 
neutrality and p e = n t kT Q we can derive that 

p e = Apl with A = ^pI~ Y T q (2) 

a constant that is evaluated using reference values of p\ and 
T & (solar wind values). Note that y = 1 corresponds to 
assuming that T & is constant, and y = 0 gives a constant p e . 


Finally, Faraday’s law is used to advance the magnetic 
field in time, 


We use a cell-centered representation of the magnetic 
field on a uniform grid. All spatial derivatives are dis¬ 
cretized using standard second order finite difference sten¬ 
cils. Time advancement is done by a predictor-corrector 
leapfrog method with subcycling of the field update, de¬ 
noted cyclic leapfrog (CL) by Matthews (1994). An advan¬ 
tage of the discretization is that the divergence of the mag¬ 
netic field is zero, down to round off errors. Also, the dis¬ 
cretization conserves energy very well (Holmstrom, 2011). 
The ion macroparticles (each representing a large number of 
real particles) are deposited on the grid by a cloud-in-cell 
method (linear weighting), and interpolation of the fields 
to the particle positions are done by the corresponding lin¬ 
ear interpolation. Initial particle positions are drawn from 
a uniform distribution, and initial particle velocities from 
a drifting Maxwellian distribution. Further details on the 
hybrid model can be found in Holmstrom (2010). 

The implementation of the algorithm was done in the 
FLASH software framework, developed at the University 
of Chicago (Fryxell et al ., 2000), that implements a block- 
structured adaptive (or uniform) Cartesian grid and is par¬ 
allelized using the Message-Passing Interface (MPI) library 
for communication. Further details on the FLASH hybrid 
solver can be found in Holmstrom (2011). The general 
solver will be included in future releases of FLASH. 

A crucial point in modeling the interaction between the 
Moon and the solar wind is how to choose the inner bound¬ 
ary conditions, at the surface of the moon. Here we model 
the Moon as an absorber of the solar wind. Ideally we 
should then remove all ions that hit the surface of the Moon, 
but that would introduce a discontinuity in the charge den¬ 
sity, from solar wind values outside the sphere, to zero in¬ 
side the sphere. However, numerical differentiation of dis¬ 
continuous functions on a grid will lead to oscillations and 
instabilities in the fields. Therefore we choose to keep also 
ions inside the Moon, but we gradually reduce their weight 
by a factor f 0 b S after each time step while they are inside the 
Moon (the mass and charge of each particle are reduced). 
Thus, the Moon acts as a gradual sink for ions, and this ap¬ 
proach also avoids the problem of pi = 0 in Eq. (1). How¬ 
ever, we do not plot these reduced ions in what follows, only 
protons that have never hit the Moon are used when com¬ 
puting the various statistics. This can be seen as a form of 
smoothing, done in a self consistent way. Often smoothing 
of fields are used in hybrid solvers, but it is then done in an 
ad-hoc way, as an extra step. The downside of our approach 
is that we introduce a numerical parameter, but experience 
shows that the solution will not change by much when re¬ 
ducing / 0 bs, indicating that the solution is converging. The 
reduction parameter cannot be chosen too small however, 
since a steep gradient in ion density will cause oscillations 
in the magnetic field. We have verified numerically that this 
approach keeps the IMF constant inside the Moon, and thus 
that no spurious currents are present there. There are no 
inner boundary conditions on the electro-magnetic fields. 

The coordinate system is centered at the Moon, with the 
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Fig. 1. The proton number density for an IMF along the direction (1, 1,0). This also illustrates the coordinate system and the simulation geometry. The 
cuts are the planes x = —10000 km, x = —30000 km, y = 0, and z — 0. These cutting planes are the same that later are used in Figs. 2-4. 


a- axis is directed toward the Sun, so the solar wind flows 
opposite to the a- axis. The IMF is in the xy-plane, and 
we will from now on denote it as the IMF plane. This 
defines a plane of symmetry, since the only asymmetry in 
the problem stems from a non-zero component of the IMF 
perpendicular to the x-axis. The solar wind consists only 
of protons (and the massless electron fluid), and we have a 
rectangular simulation domain divided into cells. The +x- 
face of the domain is an inflow boundary where we have a 
layer of shadow cells outside the domain that are filled with 
macroparticles drawn from the solar wind proton velocity 
distribution, and where the solar wind magnetic field is 
set. The —x-face is an outflow boundary where particles 
exiting the domain are removed, and where the magnetic 
field is extrapolated from the interior. The four other faces 
have periodic boundary conditions for fields and particles. 
The simulation domain is divided into a Cartesian grid with 
cubic cells of size 160 km. The geometry of the simulation 
domain is shown in Fig. 1. The extent of the simulation 
box is [-32.40,4.08] x [-8.32, 8.32] x [-8.32, 8.32] • 10 3 
km ~ [—18,2.4] x [— 4.8, 4.8] 2 7?l, with the Moon as a 
sphere of radius = 1730 km at the origin. The total 
number of simulation macroparticles is about 300 million, 
with 121 particles per cell on average. Here we have used 
values of y = 5/3, f 0 b s = 0.965, a time step of 0.08 s, 
and a final time of 170 s, when the solution has reached a 
steady state. The plots in what follows show the solution at 
the final time. Five subcycles of the CL algorithm are used 
here when updating the magnetic field, and the execution 
time is about 13 hours on 384 CPU cores. 


3. Results 

In what follows, we have used typical solar wind condi¬ 
tions at 1 AU (table 4.1, Kivelson and Russell, 1995), un¬ 
less otherwise noted. The solar wind velocity is 450 km/s, 
the number density is 7.1 cm -3 , and the ion temperature is 
1.2 • 10 5 K. The electron temperature is 1.4 • 10 5 K, and 
the IMF is in the xy-plane with a magnitude of 7 nT. For 
these plasma parameters, the ion inertial length is 85 km, 
the Alfven velocity is 58 km/s, the thermal proton gyro ra¬ 
dius is 66 km (the gyro time is 9 s), and the ion plasma beta 
is 0.6. 


3.1 Overview 

In Fig. 2 the number density (three first rows) and the 
ion velocities (three last rows) are shown, each for three 
different IMF directions (at 0, 45 and 90 degrees to the 
x-axis, in the xy-plane). The most obvious feature of the 
interaction between the Moon and the solar wind is that a 
wake is formed behind the Moon, visible as a region of low 
density. We also see how the refill of the wake depends 
on the IMF direction. The wake refills along the magnetic 
field lines, making the wake shorter when the component 
of the IMF that is perpendicular to the solar wind is larger 
(Fig. 2(1)). When the IMF is anti-parallel to the solar wind 
flow there is hardly any refill of the wake over the length 
of the simulation domain, also the flow then is cylindrical 
symmetric around the Moon-Sun axis. That the wake refill 
occur along the magnetic field lines is also visible in the 
cuts of the far wake (Fig. 2(ld)), with the refill region as a 
rectangle aligned with the IMF. 

Another notable feature is the density rarefaction cone 
that emanate from the terminator region. This region of 
low density propagate mostly above the poles in the case 
of a perpendicular IMF (Fig. 2(1)), indicating that it prop¬ 
agates perpendicular to the magnetic field lines. This will 
be discussed in more detail when we look at the associated 
reduction in magnetic field strength (Fig. 3). 

We now examine the ion velocities, computed from the 
ion current and the charge density as Ji/pi, as is usual for 
particle in cell models, since the ion current density and the 
ion charge density are quantities that are available on the 
computational grid. The small white region downstream 
of the Moon is a region without any particles that has not 
hit the Moon, thus no velocities can be computed there. 
The x-velocity is shown as the deviation from the solar 
wind velocity, when the IMF is perpendicular to the flow, 
in Fig. 2(4). We see a slight slow down of the solar wind 
in the rarefaction cone. There is also a slow down behind 
the +z hemisphere (and a speed up behind the — z hemi¬ 
sphere) seen in Fig. 2(4b). This is probably a geometric 
effect, where part of the velocity distribution function is 
emptied by collisions with the lunar surface. The y-velocity 
(Fig. 2(5)) shows how the wake symmetrically refills from 
both sides of the wake. The ions move inward not only in 
the central part of the wake, but also in the rarefaction cone. 
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Fig. 2. Proton number density and velocities in different planes, for different upstream IMF conditions. Arrows show the direction of the IMF. Top 
three rows show number density for different upstream IMF directions. Along the directions (x, y ) = (0, 1), (1, 1), and (1, 0), respectively. The 
geometry is repeated in the lower three rows (4-6) for the proton velocities u x — w sw , u y , and u z . The columns are from the left, cuts in the planes 
z = 0 (seen from +z), y = 0 (seen from +y), x = —10000 km, and x = —30000 km (seen from +x). Figure 1 provides a three-dimensional 
illustration of the cutting planes. The two vertical lines in (l-3b) show the position of the cuts perpendicular to the x-axis. 


For the case of an aligned IMF, the z-velocity (Fig. 2(6)) 
also shows a flow in the rarefaction cone toward the wake 
center. There we also see a flow of ions at the wake bound¬ 
ary in the IMF plane, and there seem to be a small kinetic 
effect of alternating upward and downward going ions in the 
upper and lower part of the wake (Fig. 2(6b)). This seems 
associated with the wave pattern in density inside the rar¬ 
efaction cone (Fig. 2(1-3)), possibly caused by oscillations 
at the wake boundary. 

The magnetic field magnitude (three first rows) is shown 
for three different IMF directions in Fig. 3, where also each 
of the vector components (rows 4-9) are shown for two 
different IMF directions (perpendicular and aligned to the 
solar wind flow). 

The most prominent feature is the enhanced magnetic 
field in the wake. For a perpendicular IMF (Fig. 3(1)) this 
region spreads out in the IMF plane, in the same region 
where we in the density saw the wake refill. For the aligned 
IMF the magnetic field increase is larger, and confined to 
the narrow wake (Fig. 3(3)). There are also structured fluc¬ 


tuations in the far wake region in the IMF plane. Then there 
is a cone of lower field strength in the case of an aligned 
IMF (Fig. 3(3)), corresponding to the cone of reduced den¬ 
sity (Fig. 2(3)). For the perpendicular IMF the cone is par¬ 
tial, and only extends perpendicular to the IMF, toward +z 
and — z (Fig. 2(1)). Between the central region of enhanced 
field and the cone of lower field, there are field fluctuations 
parallel to the cone, as is also seen in the density (Fig. 2(lb, 
2b, 3b)). Examining these density plots, the waves seem to 
originate at the wake boundary and could be caused by os¬ 
cillations at this interface between high and low ion density. 

The magnetic field z-component for the case of a per¬ 
pendicular IMF (Fig. 3(8)) exhibits a pinching of the field, 
consistent with earlier predictions (figure 3, Owen, 1996) 
and model observations (Kallio, 2005). The field is bent to¬ 
ward the IMF plane in the wake (Fig. 3(8c, 8d)). However, 
the x -component of the magnetic field is also perturbed 
(Fig. 3(4)). Close to the IMF plane, the field is bent toward 
the moon on one side of the wake, and away from the moon 
on the other side of the wake (Fig. 3(4a)). Note also that 
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Fig. 3. The magnetic field magnitude, and components in different planes, for different upstream IMF conditions. Arrows show the direction of the IMF. 
Top three rows show magnetic field magnitude for different upstream IMF directions. Then the magnetic field x-component for an IMF perpendicular 
to and opposite the solar wind flow is shown on row 4 and 5. Similarly, the y-component is shown on row 6 and 7, and the z-component on row 8 
and 9. The geometry of the cuts in the different columns are the same as in Fig. 2. The colorbar is for the row immediately to the left, and for all 
following rows, until the next colorbar. 


these bends in the magnetic field are not confined to central 
parts of the wake, but expand along with the rarefaction in 
magnetic field (Fig. 3(ld)) and density (Fig. 2(ld)). Thus, 
for a perpendicular IMF the magnetic field in the wake is 
not only bent toward the IMF plane, but also toward the 


moon, as shown in Fig. 4. 

The Mach cone of decreased density and magnetic field 
makes an angle to the solar wind flow direction of about 
10 degrees, as seen in Fig. 2 and Fig. 3. Since the outer 
edge of the cone is diffuse, it is impossible to give a 
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Fig. 4. Magnetic field lines 10000 km behind the Moon for an IMF along the y-axis (same as in Fig. 3(1)). (a) seen along the +x-axis, and (b) along 
the —z-axis. The color scale shows the magnitude of the magnetic held similarly to Fig. 3. To enhance the bending of the held the magnetic held 
x-component was multiplied by 4, and the z-component by 2, before drawing the held lines. 
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Fig. 5. The electric held magnitude, and the three components in different planes, for an IMF that is perpendicular to the solar wind how direction. 
Arrows show the direction of the IMF. For the magnitude and z-component, the solar wind convective electric held, E sw = —u sw x B sw , has been 
subtracted. The cuts in the different columns are the same as in Figs. 2 and 3. 


more exact angle. The ion-acoustic wave velocity is v$ ~ 
^/Tyk^TJnii ~ 58 km/s, and the Alfven velocity is also 
la ~ 58 km/s for these solar wind conditions, correspond¬ 
ing to angles of 7 degrees. We see that the density deple¬ 
tion travel perpendicular to the magnetic field. This together 
with the fact that the propagation velocity is larger than both 
the ion-acoustic and the Alfven velocity suggests that the 

fast magnetosonic wave, with a phase velocity of Jv\ + 
leading to an angle of 10 degrees, is responsible for the cre¬ 
ation of the Mach cone in the hybrid model. This is con¬ 
sistent with the observation of Wiehle et al. (2011) that all 
three MHD modes (fast, Alfven and slow) will be mani¬ 
fested in the lunar wake. 

In Fig. 5 the electric field magnitude (top row) is shown 
along with the three vector components (three lower rows) 
for the case of an IMF perpendicular to the solar wind 


flow direction. What is shown is the electric field with the 
solar wind convective electric field, E sw = — u sw x B sw , 
(directed along the z-axis) subtracted. We do not show the 
electric field for the case of an IMF aligned with the v-axis, 
since it is indistinguishable from zero using the same color 
scale. We can note that the regions of enhanced electric field 
magnitude coincide with those of enhanced magnetic field 
magnitude shown in Fig. 3(1). Also, there are similarities 
in the components of the electric and magnetic fields. E y 
is very similar to B z (Fig. 3(8)), and correspondingly E x is 
similar to B y (Fig. 3(8)). 

3.2 Comparison with WIND observations 

We here compare the hybrid model results with obser¬ 
vations by the WIND spacecraft on December 27, 1994, 
when it traversed the lunar wake approximately perpen¬ 
dicular to the solar wind flow, in the plane of the IMF, at 
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Fig. 6. The hybrid model magnetic field magnitude [nT] along y at x = —6.5 R l = —11245 km for two different upstream solar wind conditions 
(labeled before and after), compared to WIND observations from Ogilvie et al. (figure 1, 1996). The vertical lines show the location of the optical 
shadow, and the x-axis shows time on December 27, 1994, as in the original plot. 


a distance of x = —6.57? L ~ —11300 km. Ogilvie et 
al. (figure 1, 1996) plots several plasma parameters during 
the wake crossing. Apparently the solar wind conditions 
changed during the observation. This is most apparent in 
the plasma number density that decreased by more than a 
factor of two from the undisturbed solar wind on the in¬ 
bound trajectory to the solar wind on the outbound trajec¬ 
tory. Therefore we make two runs of the hybrid model for 
the two sets of solar wind conditions, here denoted before 
and after. The solar wind conditions used for the before 
(after) case are as follows. A solar wind velocity of 470 
(500) km/s, a number density of 5 (2) cm -3 , an ion temper¬ 
ature of 0.9 (0.5) • 10 5 K, and an electron temperature of 1.5 
(2) • 10 5 K. The IMF was (6.39, 2.33, 0) nT for the before 
case, and (5.57, 5.02, 0) nT for the after case. The time step 
was 0.08 (0.05) s and f 0 b S = 0.98(0.985). Since the plasma 
parameters are different in the two cases, these numerical 
parameters have different values to ensure the stability of 
the solutions. 

In Fig. 6 we show the model magnetic field magnitudes, 
compared with the observation by WIND. The observed 
magnetic field magnitude matches well to the model values, 
if we assume that there was a change in the upstream solar 
wind conditions around wake entry (there we switch from 
comparing the before model results to comparing the after 
model results with the observation). Before the wake cross¬ 
ing there is a drop in the observed magnetic field strength 
that is also seen in the model. This is the rarefaction wave 
that we saw in Fig. 3 earlier. The IMF in the before case 
has the largest component along the x-axis and therefore 
the magnetic field magnitude should best match Fig. 3(3a), 
with a pronounced rarefaction cone. Then there is an in¬ 
crease in the central wake which is reproduced in the after 
model. Followed by a gradual decrease at wake exit, but 
no rarefaction wave. For the after case the IMF is at about 
45° away from the x-axis and thus correspond to Fig. 3(2a), 


with a much weaker rarefaction cone. In the observation 
there are some large variations further away from the wake 
that are not present in the model. This could be upstream 
solar wind disturbances. 

In Fig. 7 the proton velocity observed by WIND is com¬ 
pared to the velocity in the two hybrid model runs. 

Again we have fairly good agreement between the model 
and observations, assuming a shift in solar wind conditions 
around wake entry. There are more variations in the ob¬ 
served velocities, but the locations of the velocity changes, 
and the slopes match the model. In the central wake there 
is however no data from the model since the model den¬ 
sity is too low. This can be seen when we compare the 
proton number density from the model with observations in 
Fig. 8. We see that the general features of the observed den¬ 
sity profile across the wake is found in the model, including 
a shift of the minimum density toward the outbound part of 
the trajectory, if we assume the shift in solar wind condi¬ 
tions around wake entry. However, the drop off in density 
is quicker in the hybrid model, and is much below the ob¬ 
servations in the central part of the wake. The cause of this 
could be the incomplete treatment of electron dynamics in 
a hybrid model. Since the electrons are a massless charge 
neutralizing fluid, there are no forces from the electrons on 
the ions, except for the electron pressure gradient term in 
the definition of the electric field (Eq. (1)). In reality, the 
refill of the wake is governed by the electrons. Due to their 
higher thermal velocity they will speed ahead of the ions to 
refill the wake, setting up an ambipolar electric field due to 
the charge separation that will accelerate the ions into the 
wake (Halekas et al., 2011). 

Similarly, the lack of electron dynamics in the hybrid 
model will affect the electron temperature, T e , in the wake. 
Observations show an increased electron temperature in the 
wake (figure 1, Ogilvie et al ., 1996). But from Eq. (2), and 
the ideal gas law, T e ~ = p 2/3 for y = 5/3, so T e in 

the hybrid model is then lower in the wake. 
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Fig. 7. The hybrid model proton x-velocity magnitude [km/s] for two different upstream solar wind conditions (labeled before and after), compared to 
WIND observations. Otherwise similar to Fig. 6. The reason that the WIND velocity is double valued in the center of the wake is that the velocity of 
the two proton beams refilling the wake from each side is computed separately, as explained in the original publication (Ogilvie et al., 1996). 



Fig. 8. The hybrid model proton number density [cm 3 ] for two different upstream solar wind conditions, compared to WIND observations. Otherwise 
similar to Fig. 6. 


4. Discussion 

Since the only deviation from axial symmetry is the com¬ 
ponent of the IMF that is perpendicular to the solar wind 
flow, we can see that the two extreme cases are when the 
IMF is perpendicular to the flow, and when the IMF is 
aligned with the flow. For a constant IMF magnitude, the 
flow and field for cases when the IMF is at an angle to the 
flow will be a mixture of these extreme cases, as seen in 
Fig. 2(2) and Fig. 3(2). 

There are many processes related to electron dynamics 
that are not captured by a hybrid model. The process of 
plasma expansion into a vacuum has a lot of details, and in¬ 
stabilities that can form (Birch and Chapman, 2001). Also, 
the region near the lunar surface probably has strong electric 
fields, and charge separation effects (Farrell et al ., 2008). 
However using the hybrid approximation it is possible to 
model the global three-dimensional interaction between the 
Moon and the solar wind, so full PIC models and hybrid 


models are complementary tools in the study of this inter¬ 
action. Also, the Moon-solar wind interaction, although 
seemingly simple, is complex enough to be a good test of 
the applicability of different models; fluid, hybrid, and PIC 
models. 

We have here considered the simplest possible model of 
the Moon in the solar wind—a spherical sink of solar wind 
protons. There are many physical processes that make the 
real interaction more complex than that. Two examples 
are magnetic anomalies that deflect the solar wind and will 
cause downstream disturbances (Lue et al., 2011), and the 
reflection of solar wind protons by the dayside lunar surface 
(Holmstrom et al., 2010). 

5. Summary and Conclusions 

We have presented a fully self consistent three- 
dimensional hybrid plasma model of the interaction be¬ 
tween the Moon and the solar wind and shown the resulting 
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global ion fluxes and fields for typical solar wind conditions 
and different IMF directions. Magnetic field enhancement 
in the wake, and a rarefaction wave in density and magnetic 
field are seen in the model. The two extreme cases are that 
of an IMF perpendicular to the solar wind flow, and that 
of an IMF that is parallel to the solar wind flow. We find 
that the wake for intermediate IMF directions can be seen 
as a superposition of these two states. Some kinetic effects 
were also observed in the model. There are density and field 
variation in the far wake, maybe caused by an ion beam in¬ 
stability from the refill of the wake. Visible are also kinks 
in the magnetic field at the inner boundary of the rarefac¬ 
tion cone for an IMF perpendicular to the solar wind flow 
direction. 

Observations by the WIND spacecraft were fairly well 
reproduced by the hybrid model, showing that a hybrid 
model with the moon as a spherical sink of ions is able to 
represent many global features of the interaction. However, 
the model central wake density and electron temperature 
were different from observations, probably due to the lack 
of electron dynamics in the hybrid model. 
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